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ABSTRACT 


Alternative methods for the rotor-dynamic and sensitivity analysis of 
large rotor systems are examined. The methods are assessed for their 
ability to utilize accurate models of reduced size along with effective 
procedures for describing the dynamic behavior of the systems. 

Frequency response-based techniques are developed for determining the 
steady state response to imbalance of the SSME turbopumps and the related 
eigenvalue problem. In these techniques, the rotor and housing are 
represented by reduced receptances associated with their coupling points. 
The housing may be described by all of its normal modes within a frequency 
- range of interest. The effects of truncated higher and lower modes are 
accounted for in an approximate manner. A reduced model for the rotor can 
be formed through an exact dynamic reduction to the degrees of freedom at 
the connection points to the housing. Alternatively, the reduction can be 
accomplished using modal representation of the nonspinning rotor. The 
reduced impedences (or receptances) of the rotor and housing are assembled 
through the impedances of the coupling elements to form the turbopump 
system. The size of the resulting system impedance (or receptance) is that 
of the number of the degrees of freedom (or forces) at the coupling points. 

A procedure is described for determining the sensitivity of the 
coupling forces to changes in the coupling elements and rotor speed of the 
turbopurnp systems. In addition, an eigenvalue sensitivity analysis 
technique is adopted for application to the systems. 

Computer programs were developed for the numerical implementation of 
the impedance and eigenvalue sensitivity formulated in this study. 
Finally, recommendations are made concerning further developments and 
requirements for other types of analysis. 
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nomenclature 


List of the main symbols 
SUPERSCRIPTS 

# the conjugate of a complex quantity 
’ indicates modified from a given variable 

a bar over a quantity indicates complex quantity 
belongs to seals 

a dot over a quantity indicates differentiation with respect to time 
T the transponse of a matrix 

SUBSCRIPTS 

b combined effects at bearings 

B refers to bearings 

c refers to turbopump case (housing) 

R refers to turbopump rotor 

I refers to intermediate (coupling) components, such as seals, 
bearing, eto. 

j particular element of a vector 
Imb imbalance locations on rotor 

k kept 

s for seals 

r reduced out 

X,Y components in the X and Y directions of the inertial coordinate system 
XYZ 

x,y components of a vector in the -x and y directions of the rotor-fixed x, 
y axes 

XZ, YZ elements in the XZ and YZ planes 



Brackets 


! ] 

square or rectangular matrix 

M 

diagonal matrix 

1 1. n 

column matrix 

[ p 

inverse of a matrix 

1 1 1 

transpose of a matrix 

VARIABLES, 

PARAMETERS AND CONSTANTS 


{a | position vectors of imbalance masses on rotor 

[a] matrix of undamped normal modes 

[c] damping matrix 

. [ S C J damping ratios of free housing 

['DpJ matrix of velocity dependent coefficients for the rotor model 

N r J damping ratios of free-free rotor 

jF } force vector 

[g] flexibility matrix 

i imaginary number /-I 

[k] stiffness matrix 

[m] mass matrix 

{ q } vector of generalized coordinates 

|p D l vector of imbalance forces 

a 

{ R } physical displacements vector 

[t] transformation matrix to reduce coupling forces on housing and 

rotor in terms of those on the rotor 

[Y ] receptance (dynamic flexibility) or generalized receptance matrix 

[z] impedance (dynamic stiffness) matrix 

A complex eigenvalue of coupled system 


v 



NJ 

I 

M 

{R i> 


matrix of frequencies of undamped housing or nonapinning rc^or 

rotor spin speed about Z axis, rad/sec. 

complex eigenvector of turbopump system 

displacement vector of the coupling points on the rotor 

displacement vector of the coupling points, and of the points 

of application of the imbalance forces on the rotor 
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1 . INTRODUCTION 


1 . 1 Background 

The dynamic analysis of models of large rotor systems, such as the 
turbopumps of the SSME (Space Shuttle Main Engine), may involve excessive 
computations and large round-off errors. Use of reduced size models for 
the systems is, therefore, highly desireable. To that end, the more 
frequently used methods include modal representation [1], substructure 
techniques [23 and static reduction techniques [3]. 

For large complex rotor systems, various modeling and analysis tech- 
niques vary in their ability co accurately describe the systems' behavior. 
This ability depends mainly on the configuration of the systems analyzed 
and the forcing conditions as well as the particular results sought. For 
the analysis of the turbopumps of the SSME, certain distinct character- 
istics should be considered. Accurate modeling of the housings requires 
representation by relatively large number of degrees of freedom. Also, the 
housing is unsymmetr ical about the spin axis of the rotor. In addition, 
engine test results have tended to show the presence of significant coupled 
rotor/housing modes. The coupling occurs through nonconservative velocity 
dependent gyroscopic and viscous forces and couples. Finally, experience 
has shown that predicted onset speeds of instability tend to be remarkably 
sensitive to small changes in the modal representation of the housing. 

Clearly, a strong need exists for a more thorough and careful con- 
sideration of the representation of the housing structural dynamics model 
and of the rotor/housing coupling procedure. For stability analysis the 
ability to predict how changes in the coupling parameters effect the com- 
plex eigenvalues of the assembled rotor/housing turbopump is equally 
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Important. An understanding of the nature of ohanges in the system eigen- 
values due to ohanges in the model parameters can be helpful in quiding an 
experimental/analytioal effort toward resolving any discrepancies between 
test and analysis results. Calculated parameter sensitivities would also 
be useful In evaluating design alternatives and in handling uncertainties 
in the input data. 

1 ,2 Objectives and Scope 

This study was mainly initiated for the purpose of developing improved 
dynamic analysis and modeling of large rotor systems. Specifically, the 
main objectives of the study were to: 

(1) address the problem of structural dynamic modeling of the large 
housing of I'he SSME turbopumps and its coupling to the rotor. 

(2) develop accurately reduced models and efficient procedures for 
determining the response and stability of the turbopumps of the 
SSME. 

(3) develop and/or apply selected sensitivity analysis techniques to 
determine the effect of variations in coupling parameters on the 
stability and forced response of the turbopumps. 

The main results of the study are outlined in this report. Following 
a comprehensive assessment of modeling and analysis methods for large rotor 
systems, alternative techniques were developed for the study of the steady 
state response and stability as well as parameter sensitivity of the SSME 
turbopump systems. The analysis is- limited to being linear and the 
response considered only concerns steady state conditions under rotor 
imbalance. 
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2. REVIEW OF ROTORDYNAMIC ANALYSIS OF LARGE SYSTEMS 


2.1 Modeling and Analysis 

Various Analysts have proposed and utilized different procedures in 
obtaining the response and stability of large order rotor systems. Tho 
motion of the systems was described by physical coordinates, generalized 
coordinates or combinations of both. 

Analysis strategies may be recognized as falling under one of two 
basic classes, Those using the complete system and those using the 
individual components of the system, together with appropriate coupling 
procedures. For large systems, both strategics call for reduction of the 
size of the models involved, An assessment of existing modeling and 
coupling techniques as related to analysis procedures used is presented in 
what follows. Also, a brief account is made of the associated eigenvalue 
and response sensitivity to changes in system parameters. The review is 
restricted to linear rotor systems. 

2.1 a Complete system-based methods 

Starting with a full scale model of the total system, several methods 
have been devised for reducing the model's size for analysis purposes. A 
static reduction technique was introduced by Guyan C <4 3 in which a given 
system is described by a selected subset of "master" degrees of freedom. 
The remaining degrees of freedom, called "slaves" are eliminated. The 
elimination is achieved through a transformation matrix formed from the 
associated static problem. The masters are capable of accurate representa- 
tion of the lower modes of the system. However, the master coordinates 
must be carefully selected otherwise some of the lower eigenvalues of the 
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system will be lost. Henshell and Ong C53 proposed an automatic technique 
for choosing the master coordinates. The criterion for the choice is based 
upon the corresponding ratios of the diagonal elements of the stiffness and 
mass matrices. A slave degree of freedom corresponds to a large stiffness 
to mass ratio. This is so since in a Guyan reduction it is Implicitly 
assumed that the mass terms corresponding to slave degrees of freedom have 
negligible effect on the mode shapes. This leads to the conclusion that 
either the corresponding masses are small or that the stiffnesses are 
large. 

Rouch and Kao [3] extended the Guyan static condensation technique to 
the analysis of simple rotor-bearing systems by accounting for gyrns^opic 
matrices. The selection of master degrees of freedom was made taking 
advantage of the knowledge of characteristic beam bending modes of the 
rotor. 1 >\ rears, as concluded by the authors, that the rotor-bearing 
systems sre amenable to this type of reduction because of their ‘nominally 
one-dimensional configuration. No assessment of the resulting accuracy 
with various types of analysis was given. Nordmann [6] attempted to mini- 
mise the uncertainties in the selection of the master degrees of freedom in 
the static condensation technique. This was achieved by applying the 
reduction technique to an arbitrarily divided rotor system and then 
assemble the reduced substructures to form a reduced system. The procedure 
is very laborious and no guarantee of accuracy is apparent. 

Downs [7] proposed a reduction method which produces frequency depen- 
dent mass matrices. Strategies for selection of master coordinates were 
given which allow for progressive improvement in the selection. The method 



was used to obtain the eigenvalues of a given system up to an upper 
frequency limit. 

A different approach for the reduction was recently advanced by 
several investigators. The main characteristic of the approach is to 
perform an exaot reduction by replacing the stiffness matrices in the 
static approach by that of the dynamic stiffness, or impedance matrices. 
In an eigenvalue analysis (or steady state response analysis), the reduced 
system will be dependent on the eigenvalue (or the forcing frequency) under 
consideration. Iterative procedures are then used with this "dynamic 
condensation" method in order to calculate the eigenvalue of the system. 
No approximations are made in arriving at the reduced model. This reduc- 
tion however is achieved at the expense of requiring of more involved 
computation. Several authors [8,9] proposed calculation procedures for 
determining the eigenvalues of undamped, conservatives structural systems. 

In an attempt to alleviate the problem of frequency dependence of the 
dynamically reduced models, Fricker [10] devised a method in which the 
frequency-dependent terms are retained implicitly by using dynamic 
stiffness matrices defined at a number of fixed frequencies. The dynamic 
stiffness matrices may be condensed efficiently to a relatively small 
number of "master" coordinates using a front solution algorithm. The 
method appears of promise, although it is not apparent how to utilize it in 
conjunction with nonconservative gyroscopic systems. 

A different approach for the analysis of mechanical structures is to 
describe the structures by their modal coordinates. A significant advant- 
age of representing a structure by modal coordinates is that the repre- 
sentation forms the basis for size reduction of the structure’s model. 
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This may be achieved by retaining only a small number of modes in a speci- 
fied frequency range, usually selected from among the lower modes. However, 
the particular modes which would influence the behavior of a given system 
depend on the type and location of the external loads. No general quide- 
lines are available in that connection. 

Childs [11] used undamped normal modes of a given rotor to perform 
transient rotordynamic analysis. The modes selected are those of the lower 
modes of the complete system involved. No allowance was made for the 
truncated modes. Choy et al. [12] used complex modal analysis to determine 
the unbalance response of damped rotor systems. For simplified oases 
considered, the complex modes of the system were obtained from those of the 
undamped normal modes of the system. 

2.1 b Subsystem-based methods 

A reduced size model can be obtained for a given system by first 
performing the reduction on the subsystems and then assemble them to form a 
reduced system. 

Hou [13] devised a scheme in which the subsystems involved are repre- 
sented by truncated sets of their free-interface modes. A transformation, 
derived from the displacement compatibility at the connection points of the 
substructures, is then applied to their modal equations of motion to arrive 
at a description of the -motion of the assembled system. This and similar 
procedures are labelled modal synthesis techniques. Li and Cunter [14] 
utilized Hou type approach to calculate the vibrations of large multi- 
component rotor systems. They also conducted an evaluation of the number 
of modes required to yield acceptable accuracy for the unbalance response 
of a two-spool gas turb ',ne engine [15]. In general, depending on the type 
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of structure, a significant number of modes might be necessary to yield 
response information of acceptable accuracy at specified locations of 
interest on the structure. These locations could be those where coupling 
or interface to other structural components occurs. 

Several methods exist in the literature which offer procedures for the 
improvement of analysis accuracy for a given number of retained modes. 
Some of these techniques utilize modal coordinates obtained with points of 
interest on the structure free (free-interface methods), and provide means 
for accomodating the effects of truncated higher modes in approximate 
manners. MacNeal [16] suggested a method in which the static contribution 
of the truncated modes is incorporated. This accounts for the missing 
flexibility effects (residual flexibility) due to truncation. The method 
can lead to significant improvement in the convergence in the solution for 
the dynamical problem. Childs and Bates [17] applied this technique in 
determining the transient response of a rotor and reported significant 
increase in solution accuracy for a given number of modes. 

A method, analogous to that of MacNeal, is also outlined in the 
literature to account for the truncated lower modes of frequencies below 
the frequency range of interest. The technique is described by Klosterman 
and McClelland [18], in which an approximate representation of the lower 
mode contribution to the structure response is proposed and labelled 
"inertia restraint". Further improvement of Hou's technique was introduced 
by Rubin [19]. Rubin proposed a method in which the first order approxima- 
tion of MacNeal is utilized to construct a pseudo static response of a 
substructure, acting within a structure, which includes inertial and dissi- 
pative contributions. Noah [20] has shown that inclusion of these second- 
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order effects may drastically improve: (i) the modal representation of a 
substructure for the modal synthesis of a structure* and (ii) the damping 
synthesis for the assembled system. 

A further reduction of a modal description of a subsystem, in addition 
to that achieved by imposing lower and upper cut-off frequencies, can still 
be accomplished by selecting certain modes from among those within the 
specified frequency range. Tolani [21] proposed a criterion based on the 
strain energy of a subsystem, whereby an upper limit on strain energy is 
set and only the modes having strain energy below this value are con- 
sidered, An approach for selecting the modes of a subsystem that are 
important with respect to the motion of another given subsystem is due to 
Morosow and Abbott [22], Unfortunately, the method is only applicable to 
weakly coupled systems. 

An alternative to using free-interface modal coordinates is that using 
"fixed-interface’' methods. In this method, the structure is represented by 
a truncated set of lower modes obtained with points of interest (usually 
those of the coupling points to other subsystems) on the structure com- 
pletely fixed. The retained modes are complemented by static "constraint 
modes" obtained as influence coefficients corresponding to the points of 
interest on the structure. Glasgow and Nelson [23] extended this technique 
(introduced in [24] and [25]) to nonconservative systems and applied it to 
the dynamic substructure coupling of various rotordynamic systems. A dis- 
advantage of this type of approach is that for systems with large number of 
coupling points among the components, the approach suffers the problem of 
introducing higher frequencies resulting from excessive number of con- 
straints imposed at the coupling (or boundary) points. In a transient 
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analysis, this will necessarily result in muoh smaller time increments and 
consequently will lead to excessive computational time and larger round-off 
errors. 

Complete treatment for a generalization of free-interface methods to 
nonconservative systems, analogeous to that of the fixed-interface method, 
[ 23 ], is still lacking. Craig and Chung [26] utilized first order state 
space formulation of the equations of motion to represent a given component 
by a truncated set of its complex modes. Wu and Greif [27], in a novel 
approach, utilized two successive transformations based upon consecutive 
use of free interface modes in the existing physical coordinates. The set 
of modes for the first transformation are those of the undamped, free 
interface subsystems. The second transformation is based on the dam ped 
fixed-interface subsystem in generalized co-ordinates which freezes the 
generalized interfaces. The second set of generalized coordinates is 
truncated and is supplemented by the fixed-interface coupled coordinates. 
The system eigenvalue problem is solved, after physical interface compati- 
bility is invoked in coupling the various subsystems. Further investiga- 
tion of this technique is warranted. 

Other methods applied to nonconservative systems include those using 
receptances of the subsystems. Palazzolo et al. [28,29] used hybrid 
representation of the receptances in terms of a truncated set of complex 
modes supplemented by correction matrices involving mass, stiffness and 
damping properties of the subsystems. The method wa 3 used to determine the 
steady state response of rotor -system components as well as in efficient 
eigenvalue reanalyses for the components. Applications to the synthesis of 
certain type of rotating machinery train was also made. The general 
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applicability of the method as a component mode synthesis technique was not 
demonstrated, however. 

The Dynamic reduction method discussed earlier [9] was also utilized 
by several analysts to reduce the size of the components of m system, prior 
to their coupling, to form a reduced system. Berman [ 30 ] presented a 
substructure coupling method that make use of frequency dependent response 
of the substructures (represented by reduced impedances) to interface and 
excitation forces. Leung [31] further developed his dynamic condensation 
technique for use with substructures. His method uses physical coordinates 
of the subsystems to satisfy compatibility conditions. The subsystems are 
described by a few lowest fixed-interface modes in conjunction with the 
static constraint modes of Hurty [24], The size of the system matrices 
involved is equal to the number of degrees of freedom, or masters, for each 
subsystem. A similar approach to that of Berman [30] was presented by 
Geering [32]. Geering applied his method, based on what he termed dynamic 
elasticity transfer matrix,, to obtain the transient response to periodic 
loads and the response spectra to stationary random loads. The eigenvalue 
problem was eluded to in his paper. 

A totally fresh subsystem analysis technique was introduced by Hale 
and Meirovitch [33,343- Admissible functions were employed for the repre- 
sentation of the substructures involved. The approach results in low order 
polynomial representation that simplifies computation. The admissible 
functions are not necessarily those of the substructures' eigenvectors. 
The geometric compatibility conditions are approximately enforced by the 
method of weighted residuals. In the author's opinion this method is a 
powerful and general technique for conducting dynamic substructure 
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analysis. However, the method would not allow use of experimental data at 
the substructure level. Also, no apparent means exists for applying the 
teohnique to damped nonoonservative systems. 

2.2 Sensitivity Analysis 

An important aspect of the modal representation of a linear dynamic 
system is the inclusion of accurate estimation of those system parameters 
to which the system behavior is most sensitive. The literature offers 
several studies and suggested procedures concerning such sensitivity 
analyses. 

Fox and Kapoor [35] obtained expressions for the first derivatives of 
the eigenvalues and eigenvectors with respect to design parameters of the 
original self-adjoint systems. Plaut and Huseyin [36] derived general 
expressions for the derivatives of eigenvalues and eigenvectors in non- 
self-adjoint systems. Their results can be useful in examining the 
stability behavior of nonconservative systems. Other sensitivity analysis 
techniques have also been reported which differ in their computational 
efficiency, type of dynamic system considered, and order of the derivatives 
of the eigenparameters involved. In reference [36], the determination of 
any one of the derivatives of the eigenvectors of a nonself-adjoint system 
requires use of all of the left-hand and right-hand eigenvectors. This was 
avoided in a formulation given by Garg [ 37 ], in which the calculation of 
derivatives is reduced to solving two sets of simultaneous linear algebraic 
equations for each eigensolution of interest. Rudisill [38] developed an 
alternative procedure to that of Garg which could be extended to find any 
order of derivative of the eigenvalues and eigenvectors, provided they 
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Nelson [39] discussed the computational requirements for obtaining the 
derivatives of the eigenvectors using the two approaches of [ 36 ] and [373. 
He pointed out that the first approach, in which the eigenvector derivative 
is expressed as a sum of all the eigenvectors, although analytically simple 
becomes prohibitively expensive for large systems. On the other hand, the 
second approach requires only the specified eigenvalue and eigenvector but 
calls for the premultiplication of an (n+1) x n matrix by its transpose to 
form an nth order system of linear equations. Nelson presented an alterna- 
tive simplified procedure for the oonsitivity analysis, based on the second 
approach. 

An interesting paper, in a series of publications by Simpson [ 4 0 3 , 
discusses means for readily obtaining eigenvalue and response sensitivities 
using simple modification of the basic eigenvalue method of Kron [41], The 
method employs a subsystem approach and provides for means of studying 
transmissibility of vibration between the system components. * A recent 
study by Yoshimura [42] also addresses the need of developing efficient 
sensitivity analysis of frequency response. A method was presented for 
determining design sensitivity coefficients of receptance-frequency 
response evaluative functions. The method was applied to machine-tool 
structures. 

Several sensitivity studies related to rotordynamic applications have 
been reported recently. Lund [43] presented a method to calculate 
sensitivities of the critical speeds of a conservative rotor system to 
changes in the design. Fritzen and Nordmann [44], on the other hand, 
reported on a sensitivity method applicable to nonconservative rotor system 
(involving representation of seals, oilfilm bearing, etc.). The method is 
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based on Taylor’s expansion for complex eigenvalues in terms of system 
parameters in which linear, quadratic, or higher-order formulas are 
obtained, depending on the order of c»ie derivatives retained in the 
expansion. The method was demonstrated using several rotordynamic 
applications. 


3. DEVELOPMENT OF ANALYSIS METHODS 


An SSME turbopump is a large nonoonservative mechanical structure 
consisting of a large component (housing) which is connected to a gyro- 
scopic damped subsystem (rotor) through nonconservative coupling elements 
(bearing, seals, etc.). The motion of the main components of housing and 
rotor is highly coupled. Simulation of the system's dynamic behavior is 
sensitive to the choice of the model, its size and to the parameters of the 
coupling elements. 

For this study, it is desired to develop significantly reduced size 
models for the housing which are sufficiently accurate, along with 
effective coupling and analysis procedures. The developed procedures would 
allow determining steady state response stability and sensitivity type of 
analysis. Based on a careful examination of existing techniques (of which 
an account is included in the previous section 2), alternative modeling and 
coupling techniques suitable for application to the SSME turbopumps are 
developed and are presented in this section. 

In these methods the steady state complex response of a given turbo- 
pump to the coupling forces is determined for cases of forced response to 
imbalance. The coupling forces are those at a given constant rotor spin. 
For stability analysis, the system will be moving in one of the modes of 
the system with a given complex eigenvalue. The coupling of the subsystem 
is accomplished using one of two approaches. The first approach uses 
superposition of the complex impedances of the individual components to 
obtain the impedance of the complete system. The system's impedance can 
then be used to determine the dynamic behavior of the system. The other 
approach uses the receptances of the housing and rotor together with the 
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impedances of the coupling elements to couple the system. This method is 
termed generalized receptance method since the final system's equations are 
in terms of the coupling foroes. 

The impedance, or alternatively the receptance technique, is used to 
develop a procedure for determining the response sensitivity to changes in 
the coupling parameters. 


3.1 The Models 

3. la The rotor reduced model 

The equations of motion for the spinning rotor of a given turbopump, 
while coupled to the housing, can be written as [115], (see Figures 1 and 2 
as applied to the HPOTP) 


>M 1 
L R'' J 


r xr 




r xr 


f 

XR 




. + 

K„ 



ta 






XZR 





i 1* 

/ 

-a 

./ , 


.%R. 




. 3 yr, 


. A< YR, 


1 

o 

o 

1 


r yr 

. O 

(ma x )' 


■(ma y )' 




+ <J> , 


+ 4 >* 


0 [S Jpj 


, V 


. J xz . 


. J xz , 



r" \ 
R,, 


“ m 


'R , 


f 


YR 




YR 


YR 



+ 

K„ 



a 1 



• • 


YZR 




t t 


. S XR 


. 


. 3 xr, 

‘ 

! XR . 


0 

0 


N vJj 


XR 


YR 


+ <ji ’ 


(ma y ) 


’-(ma x ) 


• i 

+ <j) < 


: j yz - 


J XZ . 


(1b) 


ORIGINAL PAG" tf? 
OF POOR QUAM7Y 



Main Combustion 
Chamber 


Coupling 
V Forces 
\ CP 


urxsc 


Housing 


Figure 1. The SSME Engine System for the Analysis Of the HPOTP 
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Figure 2. The HPOTP rotor under coupling and imbalance forces 


A dot in equations (la) and (1b) denotes differentiation with respect 


to time, t. Also 
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x, y, z 
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» Inertial coordinates system 

- Rotating coordinates system, x and y are 
fixed to the rotor 

* Displacement vectors of rigid bodies on 
rotor 

» Small rotation vectors of rigid bodies 
about X and Y axis, respectively 

« Stiffness matrices corresponding to 
motions in the X-Z and Y-Z planes, 
respectively 

|M xr } - External forces and moments, including 
cc *; ling forces 

- Products of inertia 

• Rotor's spin angle, velocity and 
acceleration 

- [a | cos<j>t -|a } sln|t 

x y 

» {a } sin<Jt + (a | oos,j>t 
x y 

- Vectors of coordinates of the imbalance 
masses measured in the rotor fixed axes 
x, and y at the respective rigid bodies 

« imbalance masses 

« diametral moments of inertia for rigid 
bodies on rotor 
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t'JpJ - polar moments of inertia for rigid body 

about the Z-axis 

[ s M R n ] » mass matrix of rotor associated with 

displacement degrees of freedom 

A reduced size model for the rotor, appropriate for steady state and 
stability analyses, oan be obtained by representing the rotor by its dis- 
placements at the coupling and imbalance locations. To that end, two 
procedure oan be used, (i) Modal representation, or, (ii) Dynamic reduction 
of the equations of motion in physical coordinates. 

(i) Modal representation of reduced rotor impedance : 

The following modal transformations are employed in conjunction with 
the equations of motion (la) and (lb) 


Kr 



X-Z and Y- Z planes, respectively. The modal matrices are obtained by 
solving the eigenvalue problem associated with equations (la) and (1b) with 
their right hand sides equated to zero. The equations in this case would 
represent the free vibration of the undamped, nonspinning rotor. The 
vectors {q XR } and { ^yr ^ are the generalized (modal coordinates) in the X-Z 
and Y-Z planes, respectively. The eigenvectors of equations (2a) and (2b) 
are normalized with respect to the inertia matrix of the rotor, or 
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[a yzk 1« NJ 




[a X zr ] - [a yzr 




( 3 ) 


in which [' V I N J is the unit matrix. Assuming symmetry about the axis of 
rotation, Z, one can write 





[k yzr ] [a YZR ] - 


L'* A 




( 4 ) 


where [ s A XR n ] and [ S A YR J are equal diagonal matrices whose elements are 
the square of the natural frequencies of the free undamped, nonspinning 
rotor. 

For steady state motion with constant spinning speed the 
acceleration <j) vanishes, hence the last term in both of equations (la) and 
(1b) will also vanish, leading to the following equations 



(5a) 


1 9 
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and 
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t" 1 ^] r 4 YR l + [ sa yrJ ^yr! " ^ a yzr^ 


YR 


\i. 


XR 


[a Y zr^ 


o 


^ A XZR^ ^xr 


+ * 2 


(ma v ) 


-J 


YZ 


(5b) 


Let 




^XR 


*YR 


’ ^XR^ 


'f A 
XR 




R 

k i 


, l p M l 


'f ' 
YR 


M 


XR 


% J 


l p X R l - 


(ma x ) 
J 


XZ 


. {Pypj } " F 


( ma v )1 


“J 


YZ 


[J ? ] 
L P J 


0 0 


0 bfcpJj 


( 6 ) 


Equations (5a) and (5b) may then be put in a block matrix form which 
includes modal damping; the modal damping is taken as in reference [^5] 


h R t * 


HrJ 0 


0 ('SrJ 


k„l * 


h*X B .] +N y0 J 


XR' 


'^^ C YrJ ^ A YrJ 


ld R } 
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rr A 1 T 1 

[A xzr j 


V 


° ^ A XZrJ ^ J p ^ ^ a yzr^ 

[a 1 t 
L Lrt YZR J J 


p 

l iRj 


_”[ a yzr ] [ J p ] U XZR ] 0 



(7) 


in which due to rotor symmetry 


also 

and 


f c xrJ " ^ c yrJ a KJ = E R A R 0 


f S ^XR'-^ " 1-'' A vr.] ” f Ar v ] 


YR< 


R S > J 


[ N E Rs J is the matrix of modal damping ratios 


(8) 


If the case under consideration is that of a steady state response to 
imbalance at a rotor speed of 5> , the imbalance forces can be written in 
terms of forward and backward whirling as 

{p R (t)| = |p R } e 1 ^ + (P*| e” 1 ^; i - Al (9) 

where {p r } is the conjugate of the vector of complex amplitudes { P R } . This 
is so since jp R } is a vector of real variables. As an example, the 
imbalance forces, associated with the HPOTP of the SSME [46] are given by 

|f Xr } - <i> 2 {ma x J cos <Jt , since {ma y } = {o } 

" f 5 }2 ma x} eUt + {5 * 2 ma x} e ’ 1,ft 

’ FxrI ^ * Pxr) ^ <>°> 
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Similarly 
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{p yR } - l 2 (maj sin Jt 


{’? * 2 ma x) eUt * {I }2 ma x) e ' i * t 


■ IV * ut + i4i a ' ut 


( 11 ) 


It follows that the coupling forces and generalized coordinates in equation 
(7) can be put in the form 


|F R 1 e 1 * 6 ♦ |F R ) e’ 1 * 6 

(12) 

fi„l e Lit * {q r } 

(13) 


Equation (7) can then be written as coefficients of and e or 


> l}t : [^I R J * 4 [D R ] - [flj K r } - [A R f |f r t P 


R J 
( 1 **a) 


where 


- lit - [-J 2 hi R j ♦ 4 [d r ] ♦ ty] r r i - [a r ] t if R + ? R ] 

(14b) 


U R ] 




f A XZR ] 


[A yzR lj 




t A YZR^ ^ J P^ ^ A XZR^ 


[ A XZR ] T tO [h v™ ] 


YZR 


NrJ 
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[°r] ■ 

_-i[-c R J [-a r J _ 

But |R r 1 - [Aj |q R l 


( 15 ) 


where {R r } are the displacements of the coupling points and of the locations 
at which the imbalance forces are applied 

(RpJ e 1 ^ - Iff"} e Ait - [A„] (|q R ) e i}t ♦ |q*} 

hence 


{R r } - [A R ] {q R l 


(16a) 



[A R ] {q* R l 


(16b) 


Equations (l^a) and (16a) yield the reduced model 

{r r 1 = [a r ] [" i R J + [d r 3 + [fi R l] [a r ] {f r + p r 1 

or 

+ (17) 

Similarly, equations < 1 4b) and (16b) lead to the conjugate representation 

|H*'} = [Y*] |F^P‘1 (18) 

However, only the reduced model of equation (17) need be determined. The 
conjugate representation given by equation (18) can then be obtained 
directly . 


23 


ORIGINAL PACK fg| 
OF POOR QUALITY 


In equation (17), the reduced impedance matrix of the rotor associated 
with the forward whirl speed (+<f>) is given by 


-1 


-1 


Kl ■ I y r] 


-1 


[lA r ] + i«D R ] + [fl R ]] U R f ( 19 ) 


in which [y r ] is the reduced receptance matrix of the rotor. 

(ii) Direct dynamic reduction ; 

An alternative approach can be utilized to construct the reduced 
impedance matrix for the rotor. This can be achieved by applying Gaui 
elimination only to the upper partition of the complete impedance matrix of 
the rotor corresponding to the unwanted degrees of freedom. The "kept", 
coordinates will be those of the coupling points and those at which 
imbalance forces are applied. To that end, equations (la) and (1b) are 
arranged as 


KJ 


R Rr 


r • 1 
R Rr 


R Rr 

®Rr 


^Rr 

. ♦ [<] 

S Rr 

, k rr. 


• 1 

W 


, R Rk. 


'O 


0 


0 

0 

+ , 

0 

. = . 

0 

, f r. 


- P R. 


F 


( 20 ) 


where the subscripts "r 1 ' and "k ,f denote reduced out and kept coordinates, 
respectively. 

At a given rotor spin speed the steady state variables may be 
written, as before, in the form 


Fj = {F} e i<j>t + {F } e' 


Ip l ... IF* 1 
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V 


R Rr 


-K 



R Rr 

B Rr 

83 < 

"^Rr 

i 4>t . 
e + 

_# 

B Rr 

. R Rk. 


— i 


=■*' 


. R Rk , 


. R Rk . 


-Ht 


( 22 ) 


The coefficients of the forward whirl e +i * t lead to 


where 


[-i 2 M rr + i}c rr * K pr ] [-i 2 M pl< ♦ X}0 rk * K r|< ] 

LI'* 2 M kr * ik k r * K kr) t"* 2 M kk + iiC kk * K kklj 


R r 

l R RkJ 
( 23 ) 


e — n 




Rr 


Rr J 


For convenience, equation (23) is written as 


Z Z , 
rr rk 


Z kr Z kk, 


"1 


1 S3 ' 

’ 0 ‘ 

' 


.c. 


F 

• • 


( 211 ) 


applying Gauss elimination procedure only to the upper partition of the 
impedance matrix, one gets [*17] 


r » » -i 

r 7 

J rr rk 


J kk 



' s„ 


’o 


Rr 





. - , 



R Rk. 


F 


( 25 ) 
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The rotor can now be expressed in terras of the reduced impedance matrix 
(reduced to the kept coordinates) as follows 

u’ k ] (4 1 - If! (26) 

where [z^] is the reduced imp e danc e of the rotor at a given rotor speed 
associated with the coordinates and forces of coupling to the housing and 
of imbalance. Equation (26) can be used for coupling the rotor to the 
housing model. 

3.1b The casing (housing) reduced model 

The equations of motion of the SSME housing while connected to the 
rotor, in terms of the displacements { R ss } of selected points, can be 
written as 


M SS^J ^ R ss ■ + ^ C SS^ + ^- K SS^ ^ R SS^ “ 

where the subscripts "SS" refer to the space shuttle main engine model in 
absence of the rotor of the turbopump under consideration but involving the 
nonspinning rotor Of the other turbopump (see Figure 1). The forces {f } 
are those of the unknown interaction or coupling forces with the excluded 
rotor. In other words, the missing rotor is represented by the coupling 
forces yet to be determined. 

A modal matrix of selected eigenvectors of the undamped, free 
vibration of the SSME (with [c gs ] » {F ss } set fco zero in ec 1 uat;Lon (27)) is 
used to represent the displacements in terms of the associated truncated 
set of generalized coordinates { q Q } 

{R ss ( t )} = f A sS^ f q c (t) J 
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where 


[ A SS ] ^ M SS~J ^ A SS^ " 


(29) 


use of the relations (28) and (29) with equation (27) yields 


N 0 J f%l + [ ' 2 °Q *C ,/2 J * KJ Kl 

' t A ssl T < F ss> ■ K? t p c f (30) 

where [a ] are sub-eigenvectors corresponding to the coupling points to the 
c 

rotor, and { F } are the coupling forces. 

G 

In equation (30), [''A ] is a diagonal matrix whose elements are the 

squares of the undamped natural frequencies of the 3SME associated with the 
retained modes. The matrix [ N C ] is the diagonal matrix of modal damping 
ratios. 

To obtain a reduced representation of a turbopump casing, only the 
coordinates { R q } at the coupling locations to the rotor are considered. In 
addition, if the retained generalized coordinates correspond to a range of 
excitation frequencies, the truncated higher and lower modes can be 
approximated using the associated residual flexibility [16] and inertia 
restraint [18] terms, respectively. The motion of the casing may now be 
described as 


( R c (t >! ■ [“inertJ f F c (t) l * t A o] [ q o (t)t + [° 


resid. 


] |F c (t)} 


inertia 

restraint 


modal residual 

representation flexibility 


and is associated with equation (30). 


(3D 
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It can be shown that the residual flexibility matrix [g . , 1 can be 

res Ida J 

expressed by [20] 


^resid. ^ 


UJ KJ " 1 

c du 0 du 



(32) 


in which "du" stands for "deleted upper" modes. Also, the inertia 
restraint matrix can be shown to take the form 


[°inert 1 - 72 K l 
X dL 



( 33 ) 


where "dL" denotes "deleted lower" modes, and X is the unknown eigenvalues 

of the turbopump system (coupled housing/rotor) or, in case of forced 

2 *2 

(imbalance) response, X ■ -<j> . 

For steady state under imbalance forces on the rotor at the constant 
speed <j> * the motion and coupling forces take the form 

lR 0 (t>) - |R 0 ) e 1 * 1 * |R*1 e" 1 * 1 ! ORa) 

• (F 0 (t)( - |F 0 ) e l * k + IF* } e' 1 * 6 <34b) 

M ft 

in which the complex amplitudes {r } and (r } are conjugates. Similarly 

0 c 


lq c (t)} « 1 q c } 



( 35 ) 


This leads, using equation (30), to 


fq 0 l, {q*l - ['i 2 i iif-2 c 0 a 0 1/2 J 



[A c ] T ( iF 0 | . ( 



( 36 ) 
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Similarly, using equation (3*0 > equation (31) yields 


IV ■ 

J°lnertJ + V + t°resldj] V 

(37a) 

and 



iO- 

J 0 inerJ + h’*l * [ W ]] I?*). 

(37b) 

in; which [V * ] is 

the conjugate of [y^]. 



For the purpose of determining the steady state response to imbalance 
of the rotor/housing system, equation (37a) is written as 


nu - w v . 


(38a) 


[y-[° 


inert. 


H* 0 ] 


: i 8 M 0 J + i*ra c 0 


resid, 

(38b) 


] 


where [y ] is the reduced receptance of the housing (case) associated with 

V 

the coupling points to the rotor. The vectors {r land { F } are the complex 
amplitudes of the physical coordinates and forces at the coupling points, 
respectively, associated with the e + ^^ part of the solution. No computa- 
tion involving equation (37b) is necessary since it would lead to conjugate 
amplitudes associated with e v . This amplitudes can be constructed from 
the first part, equation (37a) without resorting to extra calculations. 

The reduced impedance of the housing is given by 

U 0 ] - [*„]"’ (39) 

so that 

(90) 


V - [Z 0 ] |R 0 } 
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3.I0 The ooupllng elements ; 

The rotor and housing are interacting through bearings, seals and 
impeller reaction forces. These coupling forces can be related to the 
displacements and velocities at the connection points to the rotor and 
housing as follows (see Figure 3 ). 

Bearing and local stiffness at the housing; 

In the X~Z plane: 







(ill) 


Similar equation may be written for the Y-Z plane, so that the 
combined relations for the X-Z and Y-Z planes are 


f F 

XbC 


' R Xc 


v 

F 

Ybo 

■ 

R Yc 

■ * Ki ■ 

*Y0 



• 

R XbR 

r xr 

r xr 

F 

P Yb'R 


. r yr. 


. r yr. 


Seal and impeller Forces; 

The motion and resulting forces in the X-Z, Y-Z planes are coupled. 
For a given seal or impeller, j, the forces due to displacements of the 
housing and rotor can be expressed as 
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or 


F jXsc 

m 

i 

O. 

• 

k j 


R jXo “ R jXR 

, F jYsc . 


-" k J 



R - R 

["jYo "jYR, 


' » » 

R JXC “ R jXR 

t » 

R j Yo ~ R jYR 
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(43) 




■ R jXo' 


* R jXc 

F jXsc 


R JXo 

* [ C J 3 ] ■ 

R jYc 

F jYsc 

» 

■ C V 

r jxr 

• 

R jXR 



, r jxr 


. R JYR. 


Similar relations as those of equation (44) can be written for the forces 
at the rotor. Combining the forces on the case and housing due to all 
seals and impeller forces, one can express the resulting relations in a 
matrix form as follows. 


F X30' 


R Xc 


■ R Xo) 

R Yso 

■ - to 

R Yc 

. * [<] 

R Xo 

f xsr 

r xr 

r xr 

F YsR. 


r yr- 


• r yr. 


(45) 


Combining the relations for the bearings and seals, yields 


Ic 


IR 


[Kj] 


e \ 

R Ic 


IR 


♦ [Cj] 


‘Ic 


IR 


(46) 
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where [Kj] and [Cj] are the stiffness and damping matrices for the coupling 
or» intermediate components, i.e. the bearings and seals. 

For the steady state response to imbalance, at the spinning speed j>, 
the solution can be written as 


' R Io' 

M 

o 

I — i 



e 1 * 1 ♦ 

B IO 

R 


_ 



[ n IR J 


R r n 
^ IR J 


lRirJ 



(47) 


A similar form can be written for the coupling forces, leading to the 

r 

following relations for the coefficients of e i< ^ )t 



r 

F Ic 


• 

- ■ 


S3 

[Kj] * ii [Oj] 


F IR. 






n 





(48) 


where [z^] is the impedance of the intermediate components for a given <j> . 

— i At i- * i 

The impedance matrix corresponding to e , is the conjugate of 



3 . 2 The Analysis ; 

The system model may now be assembled from the individual representa- 
tion of its components. The assembled model can be used for a given $ to 
determine the steady state response to imbalance or the eigen-parameters of 
the turbopump systems. 




m 
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Two alternative procedures may be used to form the system's equations. 
The first approach is to (i) assemble the system impedance matrix from the 
individual impedances of the components. This can be identified as an 
impedance (dynamic stiffness) method. Another approach is to (ii) utilize 
the receptances (dynamic flexibilities) of the rotor and case together with 
the impedance of the coupling components to construct a reduced generalized 
receptance matrix for the coupled system. These two variations are out- 
lined in what follows. 

(i) impedanoe method ; 

For coupling the system components, equations ( 17 ), (40) and (-48 ) are 
used to form the system equations as 


tz; ] ♦ [z c ] 
1 11 c 

[z; i 

21 


[z! ] 

1 1 2 

[Zj 3 + [z r 3 

Z 22 R 


' 

\ 

R 


- * 
0 


c 




R_+ R t , 

R Imb 

k ^ 




(49) 


where the submatrices [z^Jare partitions of the impedance matrix [z^] of 
the coupling components, [z ] = [z T ] and {pj are the imbalance forces 


12 

acting on the rotor. The Submatrices [z^] would be inflated with columns 

and rows of zero elements if necessary to make them consistent with {r_ } + 

n 

{Rjmbl include displacements at imbalance forces, {Rj m p}- Equation 

(49) relates the complex amplitudes of the displacements at the coupling 
points of the base and rotor, associated with e*^ part of the solution, to 


those of the imbalance forces. The conjugate associated with e 


-i 5>t 


does 


not need to be computed, rather can be immediately deduced from that of 


34 


1’ 


OK POOK tiu.Mtv 


equation (-49 ). The system of Equations (49) can now be solved to yield the 
steady state response. 

The resultant naxiinum bearing reaction at a given bearing j can be 
determined as follows. Express 


AX , 


R jXc “ R jXR 


[*jX0 - R JXR) ^ + 


R jXR ^ 6 


~i<j>t 


AX. + AX* e"^ 

J J 


(aXj + AXj ) cos <Jt + i(AXj 


_# V 

AXj J sin 4>t 


(50) 


since the Rys are physical displacements, the coefficients of cos J t and 
sin l t must be real. 


Let A. - (AXj + AXj ) , Bj - (AXj - AX*) (51) 

Similarly, 

AYj = 0^ cos <j>t + Dj sin <Jt (52) 

in which 

Cj - (AYj + AY*) , Dj = (AYj - AY*) (53) 

The magnitude, AU ^ , of the relative displacement between the casing 
and rotor at the jth bearing location can be shown to be 


(au.) 2 - AX" + ay; 

J J 


+ V. cos 2 <£t + H. sin 2 <j>t (54) 

J J 
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where 


l j ■ I [A 1 + b j * c j 2 + d2) 

V J ■ 5 t*j ‘ B J 2 * C I - D f) 
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in which the maximum magnitude of the relative displacement is 

■ ( l j * t v ? + H j ) ,/2 )’ /2 (55 > 

The maximum bearing reaction force of the jth bearing follows directly as 

Mmax ' ( K b ) ( 4U jUx (56 > 

J J 

where K, is the combined stiffness at the jth bearing. 

b j 

(ii) Generalized receptance method : 

The case and rotor are represented by their receptances associated 

i At 

with the forward whirl e ^ , and reduced to their degrees of freedom at 
the coupling points. 

Using equations ( 1 7 ) and (38), the reduced models are 

(“r + v 1 " [y k ] tF * + V = {I i } 157 > 

and 


I R C 1 " h 0 l |F 0 I 


(58) 



- > j 

OF POOR OOALM 

Equations (58) and (57) can be put in a block matrix form as 


»0 



0 


’ F 

c 

Ho * S T .„ 

R Imb J 


0 

[i r1. 


1 

F + P 
L R RJ 


(58) 


On the other hand, the displacement-force relations for the coupling 
components may be written in terms of their impedances, at a given rotor 
speed $ , as 


' V 


R 

Q 

— f 

-K) ■ 


F 

R J 


• W R + S lmb- 


(59) 


in which {f r } is inflated with null values in locations corresponding to 

{pj, and [z T ] is inflated for compatibility with the mobility matrix in 
R 1 

equation (58). Substitution from (58) into (59) yields 


or 


F c 



( F 
1 Q 

) 

F ' 

R 

- w 

. M 

I F * 

l R 



■° I y rI- 




o 

o 

1 


/ ! 

1 F 
I C 

NJ + [z^j 

. 0 [yr1 . 


I 

1 F ' 
\ R 


‘[r 0 ] o " 


0 ' 
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( 60 ), 


( 61 ) 


This may be expressed as 
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' F o' 

» 

' = - [E p ] • 

■ 0 ' 

F 

l rJ 


F 

R . 


( 62 ) 


The coupling forces | f c | and|F R jmay be obtained from the above equation 
(62). However, since these forces are related, a more efficient procedure 
can be adopted by which the problem is stated only in terms of, say,|F R j 
To that end, introduce the following transformation 


RJ 


[t] {f r } 


(63) 


substituting in equation (62) and premultiplying by [t] T , yield 


[t] T [Yp][T] (r R | - - [t] t [e p ] 


0 ^ 


1 P RJ 


(6ij) 


or 


M IFrI - |Pr1 


(65) 


The above generalized receptance relations (65) may now be solved to 
obtain the coupling forces at the rotor. The bearing reactions can then be 
obtained by calculating {F q } from (63), { R G } from equation (38) and { R R } 
from equation (57). The maximum resultant bearing forces follows as 
before, using equation (56). 

3 . 2b Stability analysis-system eigenvalues 

The dynamic stability of the turbopump non-conservative system can be 
determined using either the impedance or the generalized receptance 
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formulations. For a given rotor speed | , it is desired to determine the 
eigenvalues, A's, of the coupled rotor/housing system. The forms of dis- 
placements and internal coupling forces become (replacing equations (15), 
(34), (12) and (34b)) 


{R R <fc) } - 

•V e Xt 

. (R 0 <t)| 

- (R 0 ) * xt 

(66a) 

|P R (t)| - 

(F- R | 

, |F 0 (t)| 

= {F q } e* fc 

(66b) 


where above any of the variables represents complex amplitude and A is 
the complex eigenvalue whose real part determines the stability of the 
turbopump system. 

Equation(65) is best suited for the eigenvalue analysis since the 
matrices involved are of small size. To that end, with the right hand side 
of the equation equals to zero, the coefficient matrix on the left hand 
side becomes the x~matrix whose determinant vanishes for values of X equal 
to the eigenvalues of the system. The X-matrix can be written, using 
equation (65), as 


det [Y( A ) ] = 0 

in which 


(67) 


where 


[m>] - [if h p (A)] [t] 


[y»)l 




tyx)] o ‘ 

o [*„(»a 


(68) 
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In expression (68), the matrices of the right hand side are obtained from 
equations (48), (37a) and (19)* with A replacing (i|). 

Conceptually, equation ( 67 ) can be solved possibly by using iterative 
procedures or other techniques. 

Another approach applicable specifically to the analysis of the SSME 
turbopumps can be formulated as follows. The number of retained free-free 
planar normal modes of the nonspinning rotor are selected so that it is 
equal to the number of degrees of freedom in the plane of the modes at the 
coupling points. This will result in square modal matrices [a r ] in 
equation (19). The impedance matrix of the rotor corresponding to a given 
eigenvalue, A, of the rotor /housing system takes the form 


[Z R ] " f^R-Isquare 


+ X [D R ] + [fl R ] 


r A r 1 

L R J square 


(69) 


The inverses [a r ] g quare and [a r ] g quare need be calculated only once for a 
given rotor. 

Now, equation (49) is written, in absence of imbalance forces, in the 

form 


4 

1 

'[z 0 ] 0 

[r ’ 

C 


. ° ihA 

.V 



(70) 


premultiply both sides by 



0 



one obtains 



1 

0 

0 

>1 

[Zj] + 

■ 1 

0 

'"7 
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M 
■ / ■ 
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f o 

11 

O 


. 0 E^rO. 


- 0 CZr]J 


W R. 


loJ 


( 71 ) 
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Equation (71) is now more amenable to manipulation for determining the 
eigenvalues. This is the case since the eigenvalues, A, appear in the 
submatrices [Y q ] , [Zj] and [z R ] in an explioit fashion. In addition, this 
formulation makes solution methods, including those of iterative procedure, 
computationally efficient. Methods such as those of reference [48] can be 
used in this case. 

It should be noted that similar formulation to that just described can 
be utilized for determining the steady state imbalance response. This 
would be advantageous as long as a square modal matrix, [ A R ] s q Uare » is 
employed. In case of forced response, equation (71) would take the form 


o 

o 

>H 

. L 

[zj] + 

Kj 

0 " 
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S R *■ f lmb 
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0 j 



it 

\\ 


(72) 


3 . 2c Sensitivity analysis 

The sensitivity analysis in this study concerns predicting the 
variations in the turbopumps* eigenvalues and steady state imbalance 
response due to changes in the rotor/housing coupling parameters. 

The eigenvalue sensitivity analysis was developed at an earlier stage 
of this study and is based on the formulations of the equations of motion 
of Childs [45]. In his approach, Childs utilized the modal parameters of 
the rotor /housing system coupled only through the stiffness at the bearing 
to form the full dynamic models of the turbopump involved. A computer 
program was developed by Childs based on formulations of [45] for 

Cv m b 1 '* < •** *' 
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determining the stability and forced response of the HTOTP and HTFTP of the 
SSME. 

This computer program was modified to incorporate a procedure for 
determining the eigenvalue sensitivity of turbopump systems. On the other 
hand, a response sensitivity method was developed based on the impedance 
and generalized receptance techniques described in previous sections of 
this report. An eigenvalue sensitivity analysis method can also be 
developed based on the receptance techniques but was not attempted in this 
study. 

(i) Eigenvalue sensitivity : 

Formulations of reference [ -4 5 U for the equations of motion of the 
turbopump are used for the following development of eigenvalue analysis 
procedure. In this development, the notations used in conjunction with 
formulations of [45] are basically retained for ease of reference. 

The equations describing the motion of the coupled housing/rotor of a 
given turbopump in terms of their respective free, undamped modes and the 
stiffness at the bearings can be shown to take the form 





N r J + [BAE R ] T ['-K bs ][BAE R ] -[bAE r ] T N bN ][BAE o ] 


V 

w- 


. + 





• • 


~[BAE 0 ] T [''K b J[BAE R ] [ s A 0 J+[BAE 0 ft' K b J[BAE 0 ]_ 


q c. 


(73) 

As' an approximation, two uncoupled equations of the form of equation (73) 
were assumed for the X-Z and Y~Z planes. 

[ N A R J = eigenvalues of the undamped, free-free nonspinning rotor as 
in equation (8), and 
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[ s A c J - eigenvalues of the undamped free housing of the turbopump 
considered in presence of the nonspinning rotor of the other 
turbopump, as defined previously in connection with equation 
(30) 

[BAE R ] * modal sub-matrix of the free-free rotor corresponding to the 
coupling degrees of freedom in either the X-Z or the Y-Z 
plane 

[baE q ] - equivalent to [bae r ] but for the housing in either the X-Z 
or the Y-Z plane 

[ S K. 1 « total stiffness at bearings due to bearing and local 

u 0‘'» J 

stiffness at housing. 

Equation (73) for the coupled rofcor/housing motion is solved for the eigen- 
values [ s A ] and modal matrix [^ xz ] and from a similar equation in the Y-Z 
plane for [ s A^ s ] and [ Z 1 * 

Utilizing the above modal parameters, the coupled equations of free 
motion of the complete tur'bopump system including damping, seals and 
impeller forces can be shown to take the form [45]. 


["Is] 


l q Y 


[cx] + [sex] i[CM] + [s c ] 
“U[cm] t + [sl T ) [cy] + [scy] 


J 1 4 Y, 


Hyr} + [SKX] 


r-iT 


Mo] + [sj 


- 

V 


■ q Y- 


0} ( (74) 


where the various submatrices are defined in [45] and reflect the coupling 
parameters, structural damping, spin rotor speed and gyroscopic terms. 
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It is desired to determine the sensitivity of the eigenvalues l n 
associated with equation (7*1) to the various parameters of the coupling 
elements such as the bearing stiffness. To that end, the eigenvalue 
problem associated with equation (7*1) is expressed as 

[a 2 [ s M s ] + A[C] + [k]J { q } - |0} (75) 

Following the method described in reference C*1*13, the sensitivity of a 
particular eigenvalue A n to a parameter p is determined to a first order 
approximation by the eigenvalue derivative. 


X 


n,p 



The derivative was 3hown In reference [4*1] to be 


(76) 


X 


n,p 


U n t T U„ M,p + X n C, p + K,p] lh n ] 


(77) 


where and { h n } are the nth left and right eigenvectors for the complex 

eigenvalue problem normalized 3uch that 


U n l T f 2 » n M + £] (h n j - 1 


(78) 


To a first order approximation in a Taylor series, a new eigenvalue X n 
corresponding to a change of Ap in the parameter p is 


X 


n 


no 



Ap 


or 


A n = A no + AA n 


(79) 


( 60 ) 


44 


In order to assess the degree of change in X , the variables involved may 
be normalized, so that comparison can be made between the following 


Relative change Is P ■ jp (81) 

AX 

Relative change in X„ - •? — (82) 

n x 


By way of example, suppose that it is desired to determine the 
sensitivity of the eigenvalue (or stability) of the complete rotor system 
to the magnitude of the equivalent stiffness at the ith bearings. Equation 
(77) in this case yields for the nth eigenvalue 



(83) 


The derivatives of the matrices in equation (83) may be determined taking 
into consideration equation (74). The result is 


M, 


K, 


" 0 , C,. 


m 0 


(84) 


and 


^ A X J,K t 




NvJ 


Y" J ,K 


j 


In equations (84), the indicated derivatives of [' A^J and [ N Ay S ] can 
be determined from equation (73) U3ing the method developed by Fox and 
Kapoor [35] in which the eigenvalue problem involves symmetric matrices K 


and M, or 


n ■ , ' 1 r r 
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(85) 


The derivative of the nth eigenvalue, y n , is 
,J n,p ’ I*/ t*. p " «,„] lv n ) 

in the case considered here. Using equation (73), this leads to 


( 86 ) 


[m], k - 0 


B, 




lT 


Hence 


^ A xO,K._ “ ^XZ ^’K W 
b j b j 


[bae r ] 

- [bae r ] t N b J, K)j 

[bae c ] ■ 

[bae r ] 

[ BAE ol T [ SK bJ'K„ 
B j 

[bae 0 ] 


<87) 

^XZ 


(88) 


an'd 


KJ,K„ - K Z [ R ]’K„ 


YZ 


(89) 


'j 


'j 


in which the {t}’s are the eigenvectors associated with equation (73). 
Therefore 


[ n aJ, k - - [l] 1 


[ ' a xJ,k h 0 
B j 

° ^J.K. 


[h] 


(90) 
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where [l] and [h] are the matrices of left and right hand eigenvectors of 
the original problem (74). This yields the desired new eigenvalues X n ’s as 


\ 


n 


X no + A n, K c 


AK r 


(91) 


(ii) Response sensitivity ; 

The sensitivity of steady state imbalance response to the various 
coupling parameters can best be determined using the impedance or the 
generalized receptance methods outlined earlier. 

Suppose that the sensitivity of the coupling forces exerted on the 
rotor is to be determined. To that end, equation (64) is rewritten in the 
form 





0 

[* R ]j 


[t] |f r } 



(92) 


where [z^] is the influted impedance matrix of the coupling components and 
is described in conjunction with equation (49). Since [z^] is a function 
of the coupling parameters (bearings, seals and impellers force coeffi- 
cients), one expresses equation (92) as 


[u(kj)j {F R } = {p r } (93) 

in which u(k.) denotes the functional dependence of [u] on the coupling 

J 

parameters of magnitude k,. These constants, k.*, can be the stiffness of a 

J J 

bearing, the cross-coupled damping of a seal, etc. 

the derivatives of the coupling forces |f r } with respect to a given 
parameter of magnitude kj is calculated as follows 
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( 94 ) 









- [ur 1 [uj, k [ur 1 



(95) 


The derivatives of the coupling forces for a given 4> and coupling 
parameters, as given by equation (95), indicate their sensitivity to 
changes in the magnitude of the coupling parameter k^ . 

Similarly, sensitivity of the steady state imbalance response as given 
by the impedance approach, equation (49), may be determined in the same 
fashion. Sensitivity of the acceleration response at the Kipling points 
or the reaction forces at these points could be derived directly from those 
of the corresponding displacement responses. 
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4. NUMERICAL IMPLEMENTATION AND EXAMPLE ANALYSIS 


4.1 Computer programs 

Two Computer Programs were developed based on the analysis methods 
presented in section 3 of this report. These two programs, written in 
fortran 77 » are as follows: 

(i) A modified version was developed of the computer program (ESTAB2), 
1983 update, which was developed by D. Childs of Texas A&M University. The 
modified program, (EIGNSENS) was developed by U.J. Fan also of Texas A&M. 
The program uses the equations of motion developed in [45] together with 
the sensitivity analysis of reference [44] to determine the complex eigen- 
value derivatives for the complete SSME Turbopump system. The theoretical 
background was described In section 3.2c. 

(ii) A new computer program (ESTABIMP) was written, mainly by U.J. Fan, to 
calculate the maximum resultant bearing forces of the SSME turbopurnps in a 
steady state response to rotor imbalance. The program can also yield the 
maximum acceleration levels at selected locations on the housing. The 
program, which is highly efficient, is based on the impedance and gener- 
alized receptance methods presented in section 3.2 of this report. 

4.2 Example Analysis 

4.2a Imbalance response of the HPOTP 

Test runs using the newly developed (ESTABIMP) were made to obtain the 
steady state response to rotor imbalance of the HPOTP, at various rotor 
speeds, 5> • The HPOTP model used is based on data given in reference [46]. 
In the turbopump configuration used, the rotor is that termed current rotor 
in Appendix A of [46] and is based on a model by B. Rowan. The housing is 



based on a 1 982 Rocketdyne model. For both rotor and housing, a 0.558 of 
critical damping was assumed. This value replaces the diagonal elements of 
the matrices [>E RS ] and ["*0 O of equations (8) and (30), respectively. 
The damping at the bearings was taken as 3.0 lb. sec/in. The stiffness of 
all the bearings were the same and of the constant value of 5x10 5 lbs/in. 
The values used for the local stiffnesses of the housing at the bearing 
locations are given by equation (1) of [46]. the seals and impeller force 
- coefficients and imbalance distribution were also as given in reference 
[46]. The housing free-free modes are used in an approximate fashion as in 
program (ESTAB2). Both the impendance and the generalized receptance 
approaches were used in the computational procedure as follows. 

The impedance matrices for the s rotor and housing are formed using 
equations (19) and (39), respectively. However, no account was made of the 
residual flexibility and inertia restraint terms for the housing as in 
equation (38b). The assembled turbopump model as given by equation (49) 
was used, in conjunction with a modified Gauss elimination technique, to 
determine the response amplitudes of the connection points on the rotor and 
housing. These amplitudes were used with equations (50) through (56) to 
determine the resultant maximum bearing reactions at various values of the 
rotor speeds, 4> . Table 1 shows a comparison of the results obtained for 
two selected values of <j> between the impedance method using (ESTABIMP) 
program and the method of reference [45], using the (ESTAB2) program. The 
slight differences noted in the results of the two methods can be attri- 
buted to the difference in the numerical implementation of two different 
procedures. The difference may also be due to the approximations involved 
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Table 1 . Comparison of the maximum bearing reaction forces 
calculated by the impedance and reference [453 methods 



Rection forces, lbs at 
24,000 rpm 

Impedance 

Ref. [45] 

98.9 

100.3 * 

114.5 

117.3 

44. 

47.5 

3.8 

2.8 


5 



in handling the housing modes which influence the representation of the two 
models in different ways. 

Experience with the programs showed that the impedance approach is 
more efficient than that of the approach of [*153 for a given rotor speed. 
However, when the steady state response at various rotor speeds are re- 
quired, the approach of C ^ 53 appears more beneficial to use. On the other 
hand, the generalized receptance method, although not fully tested, appears 
to be more efficient than both of the impedance method and the method of 
reference [45]. The receptance method requires only one inversion (see 
equation (17)) for assembling the system's generalized receptance matrix. 
In addition, the size of the assembled matrix is considerably smaller than 
that of the impedacne,, as can be deduced from examination of equation (65). 
4 • 2b Eigenvalue sensitivity of the HPOTP 

Few test cases were run using the newly developed program (EIGNSENS). 
A case considered is that of determining the effect of changing the 
stiffness of the first bearing on the eigenvalues of the HPOTP. The 
configuration and data for the HPOTP are those discussed in section 4.2a. 
The rotor speed is taken as | = 30,000 rpm and the assumed change in the 
stiffness of the first bearing was 20$, or 

K d = 0.80 x (K_ ) = 0.8 x (5 x 10 5 ) 

B 1 B 1 o 

= 4 x 10 5 lbs/in (95) 

The bearings’ stiffnesses are assumed to be in series with the local 
stiffnesses, K »s at the housing. The equivalent stiffness at the first 

C 

bearing location is therefore 
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(96) 
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B, (K Bi + K o, >' 

1 1 1 


so that in equation (87) for one of the planes 


[ K ] K 

b 1 ’ K B 


(K b ) 

D 1 ,K 


0 0 
0 0 
0 0 


(97) 


0 


0 

0 

0 


(98) 


The derivatives of the eigenvalues for the HPOTP System considered 
were then calculated using (EIGNSENS). Examination of the results can best 
be achieved by using the sensitivity measures proposed by Fritzen and 
Nordmann [44]. The eigenvalues of the complete turbopump system is 
expressed as 


X = ct + 
n n 


ito n 


(99) 


A relative measure of the sensitivity can be constructed by defining the 
following nondimens ional ratios, using the relations given by equations 
(91 ) , so 


Aa n /a no 
4K /K 
B 1 B 1 


Re ( 


3X r 
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and 


AlD M /lil^ 
n no 

AK /k r 

B 1 B 1 


3A 


Im ( 


n 


3K r 


w 


no 


(100b) 


As discussed in reference [44], the designation of relations (100a) 

and (100b) would enable comparison of the extent by which modifications of 

several parameters may have on different eigenvalues. Large ratios of 

equations (100a) and (100b) indicate that changes in K R will have a large 

1 

effect on both the damping and frequencies of a given complex mode. 

The results of the sensitivity analysis of the HP0TP eigenvalues for 
the case considered is presented in Table 2. The results were also checked 
by calculating directly the new eigenvalues using (ESTAB2) , with the first 
bearing stiffness as given by i*/uation (95). Close agreement was observed 
between the results of the direct calculation and those using sensitivity 
derivatives. 

In Table 2, the relative change of stiffness is negative, or 




0.2 


The reversed sign in the last column, of Table 2, therefore, indicates 
that the decrease in stiffness at the first bearing results in an increase 
in the damped frequency of the system. That increase is highest at the 
13th and 14-th conjugate modes as well as in the higher complex modes. 

The relative change in a nQ is of particular significance since it 
indicates the effect of stiffness change on a possible change from stable 
to unstable condition for the system. The results of Table 2 shows the 
highest sensitivity in this regard to be associated with modes number 
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TABLE 2. Relative changes in eigenvalues due to 
variation in bearing stiffness 


Eigenvalues 


a. 




Ao V/ a no 

ak r /k r 
B 1 B 1 


A to /to 
n no 

AK r /K r 
B 1 B 1 


1,2 

- 1.3467 

+ 275.81 

3,4 

-1 .4204 

+ 284.07 

5,6 

- 2.8753 

+ 530.76 

7,8 

- 2.6884 

+ 537.60 

9,10 

- 4.6546 

+ 690.71 

11,12 

- 3.4907 

+ 700.08 

13,14 

- 146.73 

+ 1323.7 

15,16 

- 6.6063 

+ 1865.8 

17,18 

- 9.4490 

+ 1888.0 

19,20 

- 5.5610 

+ 1904.7 

21,22 

- 9.7791 

+ 1948.4 

23,24 

- 8.6708 

+ 1961 .3 

25,26 

- 11.093 

+ 2207.3 

27,28 

- 11.167 

+ 2209.4 

29,30 

- 21 .747 

+2681 .5 

31 ,32 

- 13.863 

+ 2715.1 

33,34 

- 26.206 

+ 2878.8 

35,36 

- 14.735 

+ 2943.5 

37,38 

- 16.792 

+ 3058.5 

39,40 

- 15 . 61)8 

+ 3068.0 

41 ,42 

- 60.200 

+ 3170.3 

43,44 

- 2.1325 

+ 3375.0 

45,46 

- 17.354 

+ 3408.2 

47,48 

- 32.047 

+ 3463.5 

49,50 

- 384.47 

+ 4395.3 

51 ,52 

- 184.90 

+ 5069.9 


+ 1.86 

X 

10“* 


0 



0 



0 


+ 0.94 

X 

10~ 2 

+ 2.78 

X 

10 ”* 


0 



0 


+ 1.91 

X 

10" 2 

+ 2.12 

X 

10~ 4 

- 1.43 

X 

IO"* 


0 


+ 2.39 

X 

10" 3 

+ 0.57 

X 

10~ 2 

- 1.18 

X 

10" 1 

+ 2.41 

X 

10" 3 

+ 1.59 

X 

io "* 1 


0 


+ 0.67 

X 

10 -1 

+ 1 .32 

X 

io “ 3 

- 1 .54 

X 

10 " 11 


0 


+ 2.68 

X 

10“ 2 

To . 51 

X 

io " 3 

+ 1.81 

X 

10~ 3 


0 


+ 0.58 

X 

io ~ 2 


0 


+ 2.09 

X 

10" 2 

To . 56 

X 

10" 3 

+ 0.61 

X 

10' 2 

Tl .84 

X 

IO" 4 

+ 0.89 

X 

io " 1 

T 4.52 

X 

10~ 3 

+ 1 .02 

X 

icf * 


0 


+ 1.79 

X 

io “ 2 

+ 4.91 

X 

lO^ 

+ 1.12 

X 

io “ 2 


0 


+ 1 .94 

X 

10 H 

Tl .44 

X 

10~ 2 

+ 0.63 

X 

10" 1 


0 


- 6.15 



To . 57 

X 

io ” 1 

+ 1.11 



T2.17 

X 

io -2 

- 1 .51 

X 

10" 2 

+ 1 .66 

X 

10” 1 

-1 .8 : 

K ' 

io ” 1 

Tl .14 

X 

IO -1 
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15, 16, 19, 20, 33, 3 l \ and the last eight modes. 

associated with the relative change in an a indicates 

no 

damping of the corresponding mode is reduced. It might 
modes 15 and 16 could first become unstable. 


A negative sign 
that the effective 
be concluded that 


5. CONCLUSIONS AND RECOMMENDATIONS 

Baaed on the assessment of existing analysis procedures and on results 
obtained in the course of this study, the following conclusions can be made 
in relation to the SSME turbopumps. 

5.1 Conclusions 

(1) The frequency response methods developed in this study appear to 
constitute a highly effective tool for determining the steady state linear 
response of the SSME turbopumps. The methods, as formulated with a sub- 
system approach, demonstrate the following advantages} 

(a) the dynamic model of the housing is reduced to the degrees of 
freedom of the connections points to the rotor. Since the number 
of these point is relatively small, the size of a housing model 
is drastically reduced. 

(b) in cases where the housing is represented by a truncated set of 
its free eigenvalues, the methods will facilitate incorporating 
approximate representation of the truncated lower and higher 
modes of the housing. This is accomplished without affecting the 
size of the reduced model of the housing. 

(c) a reduced model of the rotor may be obtained directly from its 
representation by physical coordinates. In that case, no approx- 
imations are made. The rotor may also be represented by all (if 
desired) of its real modes. The resulting reduced model will 
still be of the same size corresponding to the number of degrees 
of freedom at the connection points. This representation also 
has the advantage of being based on the real free-free undamped 
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modes of the nonspinning rotor, avoiding the use of complex 
modes, 

(d) the formulation adopted in this study, which use3 individual 
models of the free components to form coupled systems, allows 
performing efficient re-analyses to determine effect of ohanges 
in the coupling elements. 

(e) frequency response information for the rotor at its coupling 
points may be obtained direotly from experimental test data. 

(ii) The generalized receptance formulations can yield accurate, re- 
duced Lamda matrices for determining all of t;he eigenvalues of the coupled 
rotor-housing systems. The accuracy of the eigenvalues will depend, how- 
ever, on the accuracy of representation of the individual subsystems. If 
the modal coordinates were used in constructing the receptance matrices, 
the number of the system eigenvalues which can be determined will be equal 
to the total number of retained modes for rotor and housing less the number 
of degrees of freedom at the connection points. Procedures are yet to be 
further developed to render the method practically efficient. 

(iii) The response and eigenvalue sensitivity relations developed here 
can be used effectively in guiding directions of design changes and in 
determining effects of estimated data errors on the analysis. 

5.2 Recommendations 

Further developments can be made based on some aspects of the proce- 
dures adopted for this study and on other representations of the 
components. In particular, it is recommended to 

(i) Develop analogous procedures for transient and other analysis 
requirements. The procedures would utilize reductions of subsystems to 
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their coordinates at the connection points and allow mixed representations 
of the subsystems. 

(ii) Further develop and optimize the generalized receptance method. 
Actual test analyses using this method should be made using corrected modal 
housing model. 

(iii) pursue use of the dynamic reduction approach in developing 
accurate methods of determining eigenvalues and eigenvectors of coupled 
subsystems . 

(iv) Develop an alternative procedure to the dynamic reduction 
approach as applied to all components of a rotor system. With the SSME 
turbopumps, the procedure could utilize static reduction method for the 
rotors and dynamic reduction of the housings. The relative merits of the 
two procedures can then be tested. The latter procedure can highly facili- 
tate computation of the coupled rotor/housing eigenvalues. 

(v) Extend the free-interface modal synthesis methods for application 
to nonconservative systems. These methods can be used in reducing the 
rotor's models for use in performing transient analyses. 

(vi) To explore using Fourier and Laplace transformations in the 
analysis of reduced rotor systems. The impedance or receptance formula- 
tions might prove useful in this connection. 
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